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NATIONAL ADVISORY COMMITTER FOR AERONAUTICS 


TECHNIPAId MEI^ORANDUM 1298 


ON the PROBLEM OF GAS FLOW OVER AN INFINITE CASCADE 
USING CHAPLYGIN'S APPROXIMATION* 

By G, A* Bugaenko 

1 , Some well-known results of Chaplygin's method (reference l) 
are first presented, For the adiahatic law of the state of the gas 

y 

when p = k P , the following relations hold : 



where p is the gas pressure, p is the gas density, V is the 
modulus of velocity, Pq and Pq are the values of p and p at 

the critical point of the flow at which the velocity becomes zero, 
k is the coefficient of proportionality, y is the ratio of specific 
heats, and a and 3 are constants. 

If the angle 6 between the velocity and the x-axis and the 
magnitude T equal to V^/ 2 a are considered, then, as was shown by 
Chaplygin, the equations of gas motion assume the form 




^ = 


2T ae 


(l^TjP Srp 


2T(lrT)P+l ^ 


( 1 . 2 ) 


where 9 is the velocity potential and ''tf is the stream function. 
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Chaplygin reduced these equations to the very simple form 

M 

5cp 

^ = - i ^ ' 

Sep K SV 


(1.3) 


hy introducing the new variahle <J and the constant 
the formulas 


T 



K = l-(2g+l)‘*‘ 
(1_t)23+1 


K defined hy 


(1.4) 


Chaplygin showed that for velocities far removed from the velocity of 
sound, the magnitude K is approximately unity and equations (1.3) 
can he integrated hy assuming K equal to 1. For K = 1, these 
equations go over into the conditions of Cauchy-Eiemann; hence, 
ixi = 6 + iO will he an analytical function of the complex variahle 
f = <p + if. 


The equation for the 
the approximate treatment 


elementary vector along a streamline in 
has, as is Imown (reference 3), the form 


dz = (ae^^ + he^‘^)d'P 



1 + ^ _ 1 - (l-’loz)'^ j 

2 ^/~ 2 o ^2 ' 2 a/ 2 a -^2 / 

(1.5) 


The complex pressure 


is given hy the equation (reference 3) 


y + iX = - e"^‘^)d<P + Pod - t„o 2 )^'^^J' dz (1.6) 

2. The steady potential flow of a gas through an infinite 
cascade according to the well-known scheme of Kirchhoiff with 
separations of the Jet is next considered. The vanes of the cas- 
cade will he assumed to be plane (fig. l). 


Telocity of the gas in the flow at infinity is denoted hy 
and its angle with the x-axis hy . ^, 3 ^; velocity of the gas in 

the jet at infinity is denoted hy "^2 ®^d its angle with the 
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x-axia "by %j 2* angle 6 between the velocity vectors and the 

x-£ixis lies within the range -2rt<0^O. 

The velocity field of the flow repeats itself for each dis- 
placement by the pitch of the cascade, that is, by the vector he“^^ 

The condition of constancy of the mass flow for steady flow of 
the gas gives 

Q, = - b sin (\+ ^xi) = P^2 ''&=2 “ (2.1) 

where Pqo-|^ is the density of the gas at infinity in the flow, . p ^2 

is the density of the gas at infinity in the Jet, and n is the 
width of the gas Jet at infinity. 


By making use cf expressions (l.l) for P, equation (2.1) can be 
represented in the form 


- h 


(‘ 


V. 


°ol 

2a 


2\P 


Bin (\+ %,]^) = Voo2 ^ 



( 2 . 2 ) 


The behavior of the function f = cp + i^i^ in the z-plane of the 
gas flow is now considered. For simplicity, the function f is 
assumed equal to zero at the critical point 0 of the flow (fig. l). 
From the relation 


dcp = ” ds 
03 


Vg dS 


it follows that on moving along the streamline \|f = 0, the function 
<p varies monotonlcally from -«> at the point E (infinity in the 
flow) to +O 0 at the point C (infinity in the Jet) and passes 
through the zero value at the critical point 0 where the stream- 
line branches. The value of the potential f = ® + i\(f at the 
critical point 0’ displaced by the period he“^^ rela-^ive to 
the point 0 is found and (fig. l) 


9(0') = J' Vg ds = J' Vg ds + _F ''s + J 

CMM'O' CM MM' M'O' 


Vg ds 


where MM' is a cut parallel to the axis of the cascade. 
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The first and last integrals mutually cancel and therefore 
as MM' approaches infinity in the flow-> the following equation is. 
obtained : 

<P (O’) = h cos (K + 0tt,x) (2.3) 


Furthermore, from the relation 



where dQ is the q.uantity of gas flowing in unit time between 
infinitely near streamlines, it follows that on being displaced 
by the pitch of the cascade the function \|/ receives an increment 
equal to Q/pq. Hence, 





sin (X + 


(2.4) 


In this manner, the f-plane with double-sided cuts along the 
half-straight lines parallel to the axis of reals (fig. 2) corresponds 
to the region of the gas flow (fig. 1). All the cuts, because of the 
rule by which the cascade was constructed, are obtained from the 
Initial one (the positive cp-aiis) by simultaneous displacement along 
verticals and horizontals at distances that are multiples of Q/Pq 

and h cos (\ + ^x)^ respectively. 


Because ou = e + 10 is an ana.lytic function of f = cp + i\p, 
the problem may be solved by relating these functions with the aid 
of a parameter that varies in the upper semicircle of unit radius, 
as in the Levi-Civita method. 


By considering the recti linearity of the cuts in the f-plane, 
the analytic function f (t) is found, which brings about the con- 
forms.! tremsfomiation of the f-plane into the semicircle t. In the 
t-plane, the flow of an ideal fluid about the boundary of the semi- 
circle is constructed. For this purpose, sources and vortices of 
strengths and intensities are located, as shown in figure 3, at the 
points tjp and that are symmetrical with respect to the circle 

and at the mirror reflection of these points in the diameter, that 
is, at the points toe. and taT^. At the origin of coordinates we 
place a sink of strength 2q (and a similar sink at infinity). 
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In the constructed flov in the t-plane, the upper semicircle of 
unit radius and the diameter of the semicircle are, of course, stream- 
lines so that the stream function ^ maintains a constant value at 
the boundary of the upper* semicircle t;j in the f -plane, this 
boundary will correspond to straight cuts. The complex potential 
of the constructed flow will have the form . 


fCt) = j^CY + iq.) log (t - tx) + (-7 + Iq.) log 

(-7 + iq) log (t - to„) + (7 + iq) log ^t - - 2iq log tj 


+ constant 


or 


f(t) = ~ 1^(7 + iq) log ^t + i - 2M j - 

(7 - iq) log ^t + ^ - 2Mjjj + constant (2.5) 

where 7 is the intensity of the vortex and q is the strength of 
the sources, and 


1 

2 


M = ^ 1 to,+ 


1 -) 

too J 


The arbitrary constant in equation (2.5) is chosen so that 
f(t) becomes zero at a certain point t = e^^, the position of 
which will be subsequently determined. 


Because the logarithm has multiple values, the upper semicircle 
of the t-plane will correspond to an f-plane with an infinite number 
of straight cuts y = kq(k = 0, il, ±2,...), where cp changes from 
k7 to +o°, as easily follows from equation (2.5) by substituting 
t = e^® (the arc of the semicircle) and t = t]^ where t}^ is the 
real amount of the interval (-1, +l), the diameter of the semicircle. 

In. this manner, the function (2.5) establishes a conformal 
mapping of the upper semicircle of the t-plane on the f-plane with 
the double-sided cuts represented In figure 2. 

The function (2.5) is used in the work of N. ,1. Akhiezer 
(reference 2) where it is obtained by successive conformal mappings; 
the f-pl8uie on the half plane, the half plane on the unit circle, 
and finally the circle on the upper semicircle. 
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The point t = 0 is carried hy the transformation (2.5) into 
f - +OB so that the radii AC and BC go over into the infinite 
segments (figs. 2 and 3). The conformal property of the transfor- 
mation breaks down at the points t = - 1 , t = + 1 , and t = e^ ^ 

(the point e^*- coirresponds to the origin of the double- aided cut 
in the f-nlane). The condition df/dt = 0 for t = ei^ gives 


7 + ^ 7 - iq 

COB e-M cos e-M 


In order to obtain the elements of the motion, an expression for 
the derivative df/dt is required that is represented in the form 


M = a (t-e^C)(t-e~^g)(t-t~^) 

dt 3t (t-too) (t--tbo“^) (t-too) (t-too"^) 

The quantities 7 and q are next determined. When a point in 
the z-plane of the gas flow is displaced by the pitch of the 
cascade he“^^, the point t, corresponding to the point in the z-plane, 
goes over from one sheet of the Eiemann surface to the next, passing 
once around the point E(t=ti,),asa result of which the function 
f(t) receives an Increment 7 + iq, as follows from equation (2.5). 
Because the corresponding increments of the functions 9 Eind 
are equal to Vocq^ h cos (\ + Qx>i) and Q./pq, respectively, the 

following equations are obtained: 


7 


q = 



■^ 2 ^ h cos (\ + 60 , 3 ^) 




2a 


VqoI b sin (\ + 8002 ) 


( 2 . 8 ) 


From the expression for 7 , it is evident, among other things, 
that 7=0 corresponds to the case where the approaching flow has a 
velocity at infinity perpendicular to the axis of the cascade. 


3. The function u;(t) is next determined. The function 
w = e + io is regular within the semicircle t and has the following 
properties : 

1. At the point 0(t = e^^), the real part of the function u) 
has a discontinuity, equal to Jt, because of the branching of the 
streamline. On the arc AO the angle 6 is equal to -n and on 
the arc OB it is equal to zero. 
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2. On the real diameter of the semicircle, the function 
(jo(t) = 0 + i0 is real because its Imaginary part a = 0 on the 
free Jets where ="^oo2' 


3. At the origin of coordinates t = 0, the function <*>(t) 
is equal to because at infinity in the Jet 0=0 and 

From the preceding discussion, it follows that the function 
uj(t) admits of analytical continuation in the lower semicircle and 
may be obtained by the Schwarz formula. Thus, 


U)(t) 




e^*y + t 
t 




e^*^' + t 
” elV - t 


dcp 


i 


log 


e^^ - t 
1 - te^^ 


(3.1) 


that branch of the logarithm being chosen that is equal to i€ for 
t = 0. If the third property is used from equation (3.1) for t = 0, 
it is found that e = -0002 • The value of too ia obtained from 

equation (3.1) by making use of the value of the velocity at infinity 
in the stream 


= i log 


1 - "boa exp i6o„2 
exp i6t»2 t rri 


whence (reference 2) 




[chacci - cos (%,j^ - ^2) 
cho^j^ - cos (0^3^ + 


arg tQ = arc tg 


sha^^L sin ®oo2 

chao^i cos - cos 


(3.2) 


The pressure of the gas on a blade of the cascade is then com- 
puted. On the forward side of the plate, which is a streamline, the 
pressure is obtained by equation (1.6) where Vc=2 is the velocity 

in the stream.. 

The back side of the plate is in the gas at rest where the 
pressure is constant and equal to 


P = Pq (1 


T )P+1 

Od2'' 
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If the fact that along a streamline dcp = df is considered, a 
formula is obtained for the complex pressure in the form. 


Y + iX = 


^ J (e't“ 


,-iw 


) df 


or 


Y - iX = 


^*2 


; 


(eiw _ gto) ^ 

dt 


(3.3) 


where the integration is taken over the upper semicircle of the 
t-plane in the clockwise direction. By considering that after 
analytical continuation in the lower semicircle the function <h(t) 
assumes conjugate values at conjugate points, the following relation 
is obtained: 


Y - IX = - 


Pq 


^2 r eiuj(t) ^ 
. . dt 


|t| = 1 


(3.4) 


where the integration is taken over the entire arc of the unit 
circle in the counterclockwise direction. 

Substituting the value of df/dt from equation (2.7) in 
equation (3.4) gives 


Y - IX = 


Pq* 3. ^oo2 r gioj(t) (t-e^^ ) (t-e~^^ ) (t-t~^)dt 

,1 (t-too) (t-tcc“^) (t-tcJ (t--^"l) 


Itl = 1 


(3.5) 


The function under the integral sign in equation (3.5) has 
three poles ^ "b — "fc qcy and t = 0, all of which lie within the 

unit circle. The residues of the function at these points are, 
respectively. 


(t oo- e^^ ) (tec - e"^^) 
(toE~ t^)(toc; - too ) 


exp i%^3_ 


(tc 


3^^) (too - e"^^) 


(too “ too) (too “ too ^) 


exp it^cfci 


ift 


oog 


!-exp 
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From equation (2.6), however, it follows that 
(too- ei«^)(t^- e-^€) ^ ^ 

(t»- ?J(V- ' 214 

Hence, the residues may be represented in the form 
2 -±^ exp(ieo,i -^ooi) 

exp(iea,i +Cfooi) 

2 iq. 

- exp iSoog 

From equation (3.5), applying the theorem on residues gives 

PqV^2 

Y - iX = - — ^ — [- 2iq exp ie^g + 

(y + iq) exp(i^l - <?ooi) - (7 - i<l) exp(i0o,;L +^ooi)] 

If the real and imaginary parts are separated. 




Y [2q sin ®°°1 " ^ exp(- 0 * 3 ^) - 


X = - 


pqV, 


“>2 


(7 cos + q sin Qx,;^) exp 0^{] (3.6) 

[- 2 q cos 6^2 + ®o°l ^ ®°= 1 ^ ®xp(- Oo^i) " 

(7 sin ©cq - <1 cos e^xi) ®xpa„^ (3.7) 

The velocity on the contour of the blade is assumed to remain 
finite; the total -pressure force on the blade will then be perpen- 
dicular to the velocity and therefore X = 0, that is, 
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(7 sin ©oa + q. cos ^1) exp.(- ) - 
(7 sin - q. C03 ^3^) exp = 2q cos 


(3.8) 


Thus the pressure force of the gas on a blade of the cascade is 
hy equations (3.6) and (2.8) equal to 


Y = - 




[- - 




[^coa(A+ 0j3^^) cos 00^3^ ^ " 


2a 

Vcol 

2a 


sin(\ + ®Col) sin 0(»2 + 


2\P 


sin (\ + 00^3^) sin ©0^3^] exp(- 0„,3^) - 


f ^Ocr 

[cos(X+ 0a=i) COS 0tsi - 1 1 sin (\+ 00^3^) sin 00 ^ 1 ] exp ^ooi V 


(3.9) 


If in equation (3.9) 3 is set equal to 0 and the 

magnitude Ooo3_^ determined by equations (1.4), is correspondingly 
replaced by 


T ocl 


ODl 


|i = 1 log lai . log ^ 


2T 


Too2 


vw2 


®2 


the formula for the pressure is obtained for the case of an ideal 
fluid (reference 2). 


In order to determine the angle entering equation (3,9), 
between the axis of the cascade and the x-axis, the ratio (2.6) 
and the values of Y and q from equations (2.8) are used. Thus 


ctg(\ + 0^^) = i 1 1 - 


V 


OdI 


2\P 


2a 


2 cos € - (M + M) 
M - M 


(3.10) 


In order to compute the length of a cascade blade, equation (1.5) 
is used. Replacing dcp by df gives 



1 - te^^ ^ t - (t - e^*')(t - e~^^)(t - t~^)dt 

e^^- t . te^^ - 1 / (t - too) (t - ice) (t - tCTT^)(t - tcc"^) 
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This expression may he put in the form 

dz = iS. g. (t)dt + 2^ g2(t)dt (3.11) 

7t It 

where 

g^(t) = e^C(t-e-i€)2.(t-t~^) 

(t-ttt) (t-ta„“^) (t-to^) (t-tg,"^) 

g2(t) = e.~.^5.(tzgl^)^(t-f 

(t-tj (t-t J^) (t-tj (t-t^“^) 

The expansions of and g2(t) into the sum of simple 

fractions are of the foim 


8u(t) = 


. _1h_ . . h; (v-1,2) 


.-1 t-^-1 t 


where 


^1 = ^2= -2^ exp(i^l - 


®1 = ^2 = exp(ieooT + 0»i ) 


2iq 


oo]^ ■+ 


En = -e' 


•ie 


Cl = Ag = exp(c^3^ - i^i) 

®1 = ®2 = e3cp(- 1^3^ 


Eg = -e 


ie 


These expressions for the coefficients are obtained if the 
following relations are used: 
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exp iiJLjjoi 


e^^-toD 


(■too-too) (-tbo-toT^) 


= y+iq. 

2iq. 


exp iojgoi 


l--^e^^ 


(tore-i«^)(^-e^^) _ 

("S»-toc) (W\x=~^) 

If equation (3.1l) is integrated over the upper semicircle in 
the t-plane in a counterclockwise direction and if relation 
z(-l)-z(l) = Z is used, the following expression for the length of 
a cascade blade is obtained: 


I = 


Tt 


Ailg 


"ill! 

”i-t, 


-i-t. 


+ Bnlg + C^lg 


-1-t, 


-1 


1-t, 


-1 


-1-t 


T -1 




+ E^lg (-1) 


ab 

3t 


Aglg 


“l“'tcc -l-tgo -1-tcc 

— + Bglg --""b" + 


-1 


D2lg 


-1-^, 

1-^cr, 


L-t, 


-1 


1-t, 


1-t. 


+ Sglg (-1) 
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qa 

It 


t ”^+l 

{A^ + Cj^)lg ■ - — + + D]^)lg 


tto -1 


tgT ^+1 

toT^-i. 


{A-^ + + E^)lg(-l) 




, too"Vl 

(A -2 + C 2 )lg ^ + 

toT -1 


^ "^+1 


(B2 + pj + (A2 + B2 + E2)lg (~l) 


to= -1 


aa 

It 




(A^^ + + D-j^)lg ■— + 


ICC (At + Ct ~ B^ — D-j ) + (At + Bn + ^1^ 


qb 

It 


El 


(An + Bn + Cn + Bn )lg — + 

E<r 


ia (A-^ + C]_ - B]_ - B^^) + (C;^ + B;j_ + Eg) nlj 

The magnitudes Ej_, Rgi “■ "that enter this equation are 

shotm in figure 4. Substituting the values of the coefficients and 
making use of equation (3.S) gffes the following expression for 
the length of a blade : 


I = + ^ (a + b) cos 0oo2 185^ + 

It Eg 

s (a * I.)[(n COS 0oni - q. sin ©ooi) exp(- + 

(7 cos 0(xi + T 3in expaooi] + 

[(7 COS 0ooi - q. sin eooq) exp(- aoo]^)" 


(7 cos ©aq + q. sin ©ooi) expOoofl+ qa sin 
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The fomula for the length of a blade in the case of an ideal fluid 
(reference 2) is obtained from the preceding equation for p = 0 
if 



<1 = -V 
y = 

^ ooq 


h sin (\ + 


cos (K + 


= Ig 


Vc»l 
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